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Abstract 



The method of continuous averaging can be regarded as a combination 
of the Lie method, where a change of coordinates is constructed as a shift 
along solutions of a differential equation and the Neishtadt method, well- 
known in perturbation theory for ODE in the presence of exponentially small 
effects. This method turns out to be very effective in the analysis of one- and 
multi-frequency averaging, exponentially small separatrix splitting and in the 
problem of an inclusion of an analytic diffeomorphism into an analytic flow. 
We discuss general features of the method as well as the applications. 
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There are several problems in the perturbation theory, of real-analytic ordinary 
differential equations (ODE), where standard methods do not lead to satisfactory 
results. We mention as examples the problem of an inclusion of a diffeomorphism 
into a How in the analytic set up, and the problem of quantitative description 
of exponentially small effects in dynamical systems. In this cases one of possible 
approaches is an application of the continuous averaging method. The method 
appeared as an extension of the Neishtadt averaging procedure ■ We begin with 
the description of the method. 

Let us transform the system 



by using the change of variables z Z(z, s). Here z is a point of the manifold M, 
u is a smooth vector field on M, s is a non-negative parameter, and the change is 
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1. 



The method 



z = u(z), 



(1.1) 
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defined as a shift along solutions of the equation 1 

dZ/ds = f(Z,s), Z(z,0) = z, 0<s<S. (1.2) 

Let the change z >— ► Z transform (|1 . 1|) to the following system: 

Z = u{Z,s). (1.3) 

Differentiating (|1.3|l with respect to s, we have: 

f(Z,s)=u a (Z,s)+dfu(Z,s) or u s = [u, /]. 

Here df is the differential operator on M, corresponding to the vector field /, 
the subscript s denotes the partial derivative, and [•, •] is the vector commutator: 
Uq\ = d ul u-2 — d U2 Ui. Putting / = £u, where £ is some fixed linear operator, we 
obtain the Cauchy problem 

u s = -[£u,u], u\ s = =U. (1.4) 

We call the system (|1.4|l averaging. The equation / = £it is crucial for our method. 
The vector field / is usually constructed as a series in the small parameter and not 
as a result of an application to u of an operator £, chosen in advance. 

A nonautonomous analog of ljl.4JI can be easily constructed. If u depends 
explicitly on t then / = £u also depends on t and i|l. H should be replaced by the 
system 

u s = (£u) t -[£u,u], u\ a =o = u(z,t). (1.5) 

Properties of the averaging system can be illustrated by the following example. 
Consider the non-autonomous real-analytic system 

z = Eu(z,t), zeM. (1.6) 

Here e is a small parameter, u is 2-7r-periodic in t. Let us try to weaken the depen- 
dence of u on time by the change z i— > Z (|1.2|) with / = £u. We put 2 

£u(z,t, s) = ^ia k u k {z,s)e lk \ er fe =signfc, (1.7) 
fcez 

where u k are Fourier coefficients in the expansion u{z, t, s) = X^fcez uk ( z i s)e lkt . 
Equation (|1.5|l takes the form 

u k = -\k\u k + ie<T k [u ,u k }-2iej: i+m=km<0<l {u l ,u m l (1.8) 
u k \ s=Q = u k , keZ. 

To have an idea of properties of this system, we skip in (|1.8J) the last term. 
The equations 

u k = -|fc|u fe + iea k [u°,u k ], u k \ s=0 =u k , keZ 

1 Such method of constructing a change of variables is called the Lie method. The corresponding 
Hamiltonian version is called the Deprit-Hori method. 
2 Such an operator £ is called the Hilbert transform. 
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can be solved explicitly: 

u k = e- Ws u k og iEa * s , (1.9) 

where g s is the time-s shift z\t=o l— ► z\ t=s along solutions of the system z — u°(z). 

The complex singularities of the functions u k o g 1 * of the complex variable £ 
prevent an unbounded continuation of the solutions l|1.9|) to all the set of positive 
s. Nevertheless, the functions ljl.9JI can be made exponentially small in e since s 
can be chosen of order ~ 1/e. 

If e is not small, the operator (|1.7|l can be used to smooth out the dependence 
of u on t. Indeed, for arbitrarily small s > the Fourier coefficients u k in (|1.9fl 
decrease exponentially fast in k even if u is just continuous in t. 

Certainly, these hewristic arguments cannot be regarded as a proof of the 
fact that systems of type (jl.Sjl can be used for averaging or smoothing. Rigorous 
statements and estimates are based on the majorant method. In this paper we do 
not go into technical details, but present general ideas and applications. 

If the vector field u belongs to some subalgebra \ m the Lie algebra of vector 
fields on M, it is natural to look for the change z t— > Z from the corresponding Lie 
group of diffeomorphisms. This means that / in l|1.2|l should be taken from \- The 
same remains reasonable in the non-autonomous case. Note that the operator (|1.7|l 
is such that if u(z, t, s) <E \ f° r anv ^ the vector field £u also belongs to \ I0r anv 



2. Applications 



2.1. Fast phase averaging: one-frequency case 

Mathematical models of various physical processes use systems of ODE which 
contain an angular variable changing much faster than other variables in the system. 
Taking the fast phase as a new time, we can rewrite the equations in the form 

z = eu(z,t,e), zeM, (2.1) 

where M is the m-dimensional phase space of the system, and e is a small parameter. 
The vector field u is assumed to be smooth and to depend on time 27r-periodically. 

It is well known that by a change of the variables it is possible to weaken 
the dependence of the system l|2.1() on time. In particular, by using the standard 
averaging method, it is easy to construct a 27r-periodic in t change of the variables 
z i — ► z* such that the equations Ij2.1|l take the form 

i* = su a {z*) + e 2 u*(z*,e) + eu(z*,t, e). (2.2) 

Here the only term in the right-hand side depending explicitly on time is eu — 
0(e K ). The natural K is arbitrary and vP(z) = J Q 27r u(z, t, 0) dt. 

Now suppose that u is real-analytic in z. Poincare noted in some example that 
in this case power series in e presenting a change of variables eliminating time from 
the equations, exist but diverge: terms at e k in these series have the order kl. In a 
general situation this statement has been proved by Sauzin |21| . 
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Neishtadt [T3] noted that in this case it is possible to obtain in (12.21) 

u = 0(e~ Q/e ), a = const > (2.3) 

(e is assumed to be nonnegative) . The method Neishtadt used to prove this assertion 
is based on a large (of order 1/e) number of successive changes of variables. These 
changes weaken gradually explicit dependence of the equations on time. Ramis and 
Schafke obtained analogous results analyzing diverging series, produced by the 
standard averaging method. 

It is known also that in general a constant A > a exists such that it is impos- 
sible to construct 27r-periodic in t change z ^ z*, such that u = 0(e~ A / £ ). This 
statement follows, for example, from an estimate of the separatrix splitting rate in 
Hamiltonian systems of type (|2.1|) with one and a half degrees of freedom. 

In this section we estimate a "maximal" a for which the estimate l|2.3|) is 
possible. 

Suppose that the manifold M is real-analytic. We fix its complex neighborhood 
Mc and denote by g* the phase flow of the averaged system 3 z = u°(z). 

Let Q be a compact in M and Vq its neighborhood in Mc- Suppose that for 
any real s such that \s\ < a and for any point z G Vq the map g ls is analytic at z 
and moreover, g ls (z) G Mc- We define the set 

U Q , a = |J g ls (V Q ). 

—a<s<a 

Theorem 1 |24j . Let the positive constants a,p,£o be such that 

(1) U Q , a CM C . 

(2) The vector field u is analytic in z and C 2 -smooth in t, e on UQ. a xTx [0, £o]. 
Then for sufficiently small Eq, there exists a lit -periodic in t real-analytic in z 

map F : V Q X T x (0, e j -> M c , Q C Vq C Vq, such that 

(a) The set Vq is open in Mc, 

(b) F(z,t,e) = z* = z + 0{s), 

(c) F transforms 12.1)) into 12.2)) and the following estimate holds: 

\u(z, t,e)\< Ce~ a l\ zeV Ql t e E p , e G [0, e ). (2.4) 

Theorem 1 means in particular, that in the case when components of the field 
u are entire functions of z, the quantity a in H2.4J1 can be arbitrary positive number 
such that for all s G [—a, a] the maps z i— > g ls (z) are holomorphic at any point 
z G Q. 

Proof of Theorem 1 is based on the continuous averaging. Namely, we solve 
the Cauchy problem Q1.5[) . with ^ defined by (|1.7|l . The required change of variables 
corresponds to the value s — a/e. The averaging can be performed inside a subal- 
gebra \ m the Lie algebra of all real-analytic vector fields on M. In particular, if 
the initial vector field u is Hamiltonian then u* = u(z, t, e, a/e) is also Hamiltonian, 
and F is symplectic. 

3 It would be more correct to call it by the first approximation averaged system, written with 
respect to the fast time. 
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2.2. Averaging: mult i- frequency case 

Consider a real-analytic slow-fast system 

x = LJ+e(u(y)+u(x,y,s)), y = e(v(y) + v(x,y,e)), x G T", y G R m . (2.5) 

Average in x of u and v is assumed to be 0(e). The frequency vector u> G R" is 
constant and non-resonant. 4 

We try to weaken the dependence of the right-hand side of the system on the 
fast variables a; by a near-identity change (x mod 2ir, y) i— > (x, mod 2tt, y,). We put 

Z= (y)< 57= (o)' ™=U)' 

Note that the vector fields 57 and w commute. The system <|2.5[) takes the form 

z —ZJ + e(w + w). (2.6) 

Al. Diophantine condition. We assume that n > 2 and the frequency 
vector lo is Diophantine: there exist 70,7 > such that 

\{k,uj)\ >7o||fc|r 7 , for any k G Z" \ {0}. (2.7) 

To formulate the next assumption, we need some definition. Let g l be the 
phase flow of the system 

z = w(y). (2.8) 

For any real-analytic function f(y) with values in C n+m and the vector k G Z" 
we put /fc = /(i/)e*' fe ' x ' . The function 

gfcj = e v % (/fc off ), sec 

does not depend on x. Here is the differential of the map g*. Note that gjj, s G R 
include shifts along 5* in purely imaginary direction. We put 

E g = {x G C7(2ttZ)" : |Im^| < q, j = l,...,n}, 

V„ = {yGC™:Reye(P + z,)cR m , |Im W | < 1/, l = l,...,m}, 

where T> is a compact domain and I? + 1/ is the ^-neighborhood of T>. 
Expand the function w into the Fourier series: 

w(z,e)= fi k (y,ey {k ' x) - 

A2. Analyticity. Let the constant a > be such that for any real s G 
[—a, a], for any k G Z and for any (x, y) G S g x V v = closure(£,j) x closure(V„) the 
map (a;, y) 1— ► g w (x, y) is analytic in x, y and the map y 1— » g%w k (y) is analytic in y. 

4 The case of non-constant frequencies (cj = u>{y)) can be reduced to this one in a small neigh- 
bourhood of an unperturbed (may be, resonant) torus {y = y }. 
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Moreover, suppose that for some real p 

|g^ fe kx[o, £o )<Mll fc ir Pe " 9 " fc "' k^O, 

The function w° is of order e and we assume that g§u>°|y x j e > < s/xo- 

The constants p,, p,o, v, q, p, 7 must satisfy some conditions |18|. Here we replace 
these conditions by more restrictive, but simple ones. 

A3, p, pq, v, q, p, 7 are positive, do not depend on e, and p > 7/(7 + 1). 

Theorem 2 ^H]. Suppose that assumptions A1-A3 hold. Then there exists 
a change of variables 

z^z. = f(z, e), f : E 29/3 x V 2v/3 x [0, e ) - C n /(2nZ n ) x C m (2.9) 

such that f is analytic in z, smooth in e, f(z,e) = z + 0(e), and the system (|2.6|l 
takes the form 

z. =qj + e(W(y.) + w.(z.,e)). (2-10) 
Let w®(y m ,e) be the average in x % of w.. Then w®(y,,Q) — 0. Moreover, 

\w. - w°\ < Cperl^e-^- 1/{l+1 \ z. e E, /2 x V u/2 , (2.11) 
where C is a constant, not depending on e and p, 

q=(l+ 1 - 1 )( llo aq^ 1 /^ +1 \ (2.12) 

If the system (|2.6t is Hamiltonian with respect to a certain symplectic structure 
fl then (|2.1U|I is also £!-Hamiltonian and the change 1)2. 9|) is fi-symplectic. 

In ^Sl we a l so present another theorem which shows that the Fourier series 
w. — u>2 can be divided into 2 parts: one is small and for another we have a sort of 
control. 

Results analogous to Theorem 2 (without estimates for a) are contained in 

2.3. Exponentially small separatrix splitting 

The phenomenon of exponentially small separatrix splitting was discovered 
by Poincare Intensive quantitative studying of the problem was initiated by 
papers 0E3j (see also |2]). The method proposed by Lazutkin with collaborators 
|12l 1111 [SI H3 is based on an analysis of the separatrices in the complex domain. 
Another method was used in 0], where direct expansions of the Poincare-Melnikov 
integral in an additional parameter are analyzed. The resurgent analysis is applyed 
to these problems in |2()j . 

The main difficulty of the problem is that the traditional Poincare-Melnikov 
method can not be applied directly. Indeed, its error has the order of square of the 
perturbation. Hence, the error considerably exceeds the expected result. 
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Exponentially small separatrix splitting can be studied with the help of the 
continuous averaging method. The main idea of to reduce the rate of the perturba- 
tion to an exponentially small quantity such that its square is much smaller than 
the result. Then the Poincare-Melnikov method is applicable. 

It is natural to measure the rate of the separatrix splitting by the area ^Jo^g 
of a lobe domain, bounded by segments I S:U of the stable (s) and unstable (it) 
separatrix such that I s and I u have the same boundary points, these points are 
homoclinic, and there are no other common points of I s and I u . 

Consider the system with Hamiltonian 

H(x,y,t) = e(y 2 /2 + {I + 2B cost) cosx), (2.13) 

where x, y are canonically conjugated variables and e > is small. Theorem 1 
implies that there exists a symplectic change of coordinates x, y >— > x, y which is 

i) close to the identity, 

ii) 27r-periodic in time, 

iii) real-analytic in a complex neighborhood of the separatrices T ± of the sys- 
tem with Hamiltonian Hq = y 2 /2 + cosx 

iv) such that the new Hamiltonian function takes the form: 

H{x,y,t) = e(H (x,y) +eH 1 (x,y,s) + exp(-c/e)H 2 (x, y, t, e)). 

Here H = y 2 /2 + cosx, the constant c £ [0, 7r/2) is arbitrary, the functions Hi,H 2 
are real-analytic in x, y in the vicinity of the unperturbed separatrices of the hy- 
perbolic fixed point (x, y) — (0, 0), smooth in e > and the function H 2 is analytic 
and 27r-periodic in t. 

The ordinary Poincare-Melnikov theory applied to this system for positive c 
gives a correct asymptotics of the separatrix splitting, [25) . The following estimate 
holds: 

^„be = f^(-£)( B « B2 ' + °(i))' 

where / is an entire real-analytic function, f(0) = 2. In [251123) numerical values of 
several Tailor coefficients of / are presented. 

Analogous results for the Standard Chirikov Map and for some Hamiltonian 
systems with 2 degrees of freedom can be found in [231 1261 [T5| . 

2.4. Inclusion of a map into a flow 

In this section we consider the following problem: to present a given self-map 
of a manifold M as the time-one map (the Poincare map) in some ODE system, 
generated by a periodic in time vector field. 

The problem can be formulated for various classes of maps and the correspond- 
ing vector fields. For example, it is possible to consider generic maps and vector 
fields, reversible ones with respect to some involution, Hamiltonian, preserving a 
volume, etc. Here we discuss the analytic set up i.e., assuming that the map is real- 
analytic, wc look for its inclusion into a flow generated by a real-analytic vector 
field. The problem in C°° category is much simpler. 
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The following construction is well-known. Given a diffeomorphism T of a 
manifold M onto itself consider the direct product M x [0, 1] with the vector field 
d/dt, where t is the coordinate on [0, 1]. The map T generates the identification 

Mx{0}~Mx{l}, (2,0)~(T(z),l). 

This identification converts M x [0, 1] into a manifold M. Let 7r : M x [0, 1] — > 
be the natural projection. The smooth vector field d/dt generates on the surface 
7r(M x {0}) C M. the Poincare map coinciding with T. 

This construction does not solve the problem we deal with because in general 
it is not clear if A4 is real-analytically diffeomorphic to M x T 1 . Nevertheless, 
sometimes this can be proven. 

The problem is solved in the symplectic set up for maps which are close to 
integrable, 00121 and for generic maps |27) . 

Note that all the known proofs in the analytic set up use essentially the Grauert 
theorem on the inclusion of an analytic manifold into Euclidean space (or modified 
versions of this theorem). Our result is based on the method of continuous averaging. 

First, note that any map which is not isotopic to the identity 5 obviously can 
not be included into a flow. Let M be an m-dimensional compact real-analytic 
manifold. Let ^ be a closed subalgebra in the Lie algebra (£, [,]) of all analytic 
vector fields on M. 

We denote by X the subset of all analytic diffeomorphisms of M obtained as 
a result of the time-one shift along solutions of a system 

z = u{z,t), u(-,t)e x , te[0,i], zeM (2.14) 

(it is not assumed that u(z, 0) = u(z, 1)). We assume that the vector field u is C 2 - 
smooth with respect to time. This smoothness condition is technical. For example, 
it can be replaced by continuity in t in Hamiltonian case and in the general one. 

Obviously, all diffeomorphisms from X are isotopic to the identity inside X . 

Theorem 3 For any map T € X there exists a vector field 

U = U{z,t), U(-,t)ex, ten, zeM 

which is analytic in z and t, 2tt -periodic in t, and such that the time-2ir shift along 
its trajectories coincides with T. 

As a corollary we obtain a possibility of the inclusion of analytic maps into 
analytic flows in general, symplectic, and volume-preserving cases. 

The vector field U is obviously not unique. 

It can be also proved that if T is reversible with respect to some involution 
/ : M — > M, I 2 = idjvf (i.e., To J = I o T^ 1 ), the corresponding vector field U, can 
be also regarded /-reversible: U(z,t) = —dIU(Iz, —t). 

5 Two smooth maps Tj : M' — » M" , j = 0, 1 (M' and M" are manifolds) are called isotopic 
if there exists a family of maps T s : M' — » M" of the same smoothness class continuous in the 
parameter s £ [0, 1], such that To = To and T\ = T\. In other words, if To can be continuously 
deformed into T\. 
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Suppose that the map T is close to To (dist(T, To) = e in a complex neighbor- 
hood of M), 6 where To is included into the flow generated by a periodic analytic 
vector field U . Then the vector field U can be chosen close to U (\U — U \ = O(e) 
in a smaller complex neighborhood of M). In particular, in the symplectic case if 
T is close to an integrable map, a Hamiltonian system associated with T also can 
be chosen close to an integrable one and the orders of closeness are the same. 

Continuous averaging in the proof is used to smooth out the dependence of 
the original vector field l|2.14[l on time, |17|. 
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